WL-TR-97-3083 

FEASIBILITY  ASSESSMENT  AND 
OPTIMIZATION  STUDY  OF  SMART 
ACTUATION  SYSTEMS  FOR 
ENHANCED  AIRCRAFT  MANEUVER 
PERFORMANCE 

K.  APPA 
J.  AUSMAN 

NORTHROP  GRUMMAN  CORPORATION 
MILITARY  AIRCRAFT  SYSTEMS  DIVISION 
ONE  HORNET  WAY 
EL  SEGUNDO,  CA  90250-3277 

N.  S.  KHOT 

FLIGHT  DYNAM  ICS  DIRECTORATE 
WRIGHT  LABORATORY 
AIR  FORCE  MATERIEL  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7562 

JULY  1997 


FINAL  REPORT  FOR  JUNE  1996  TO  JUNE  1997 


FLIGHT  DYNAMICS  DIRECTORATE 
WRIGHT  LABORATORY 
AIR  FORCE  MATERIEL  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7562 


NOTICE 


WHEN  GOVERNMENT  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  ARE  USED 
FOR  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A  DEFMTE 
GOVERNMENT-RELATED  PROCUREMENT,  THE  UNITED  STATES  GOVERNMENT 
INCURS  NO  RESPONSIBILITY  OR  ANY  OBLIGATION  WHATSOEVER.  THE  FACT  THAT 
THE  GOVERNMENT  MAY  HAVE  FORMULATED  OR  IN  ANY  WAY  SUPPLIED  THE  SAID 
DRAWING,  SPECIFICATIONS,  OR  OTHER  DATA,  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION,  OR  OTHERWISE  IN  ANY  MANNER  CONSTRUED,  AS  LICENSING  THE 
HOLDER,  OR  ANY  OTHER  PERSON  OR  CORPORATION;  OR  AS  CONVEYING  ANY 
RIGHTS  OR  PERMISSION  TO  MANUFACTURE,  USE,  OR  SELL  ANY  PATENTED 
INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


THIS  REPORT  IS  RELEASEABLE  TO  THE  NATIONAL  TECHNICAL  INFORMATION 
SERVICE  (NITS).  AT  NTIS,  IT  WILL  BE  AVAILABLE  TO  THE  GENERAL  PUBLIC, 
INCLUDING  FOREIGN  NATION. 


THIS  TECHNICAL  REPORT  HAS  BEEN  REVIEWED  AND  IS  APPROVED  FOR 
PUBLICATION. 


Narendra  S.  Khot 
Project  Engineer 
Design  and  Analysis  Branch 


- - 

Mark  Hopkins,  Principal  Engineer 
Vibration  and  Aeroelasticity  Branch 


Nelson  D.  Wolf,  Chief 
Design  and  Analysis  Branch 
Structures  Division 


IF  YOUR  ADDRESS  HAS  CHANGED,  IF  YOU  WISH  TO  BE  REMOVED  FROM  OUR 
MAILING  LIST,  OR  IF  THE  ADDRESSEE  IS  NO  LONGER  EMPLOYED  BY  YOUR 
ORGANIZATION,  PLEASE  NOTIFY  WL/FIBD  BLDG  45,  2130  EIGHTH  STREET,  SUITE  1, 


DEPARTMENT  OF  THE  AIR  FORCE 

WRIGHT  LABORATORY  (AFMC) 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO 


14  OCT  97 


tV- 

X 


FROM:  WL/FIBD 

SUBJECT:  Errata  to  WL-TR-97-3083  titled,  “Feasibility  Assessment  and  Optimization 
Study  of  Smart  Actuation  Systems  for  Enhanced  Aircraft  Maneuver  Performance.” 


The  unedited  copy  of  the  subject  technical  report  was  inadvertently  distributed  in  Aug 
1997.  If  you  have  received  the  copy  of  this  TR  previously  please  discard  it. 


NARENDRA  S.  KHOT 


Research  Aerospace  Engineer 
Design  &  Analysis  Branch 
Structures  Division 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  at/erage  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  r^ucing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway.  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 


3.  REPORT  TYPE  AND  DATES  COVERED 

Final,  6-16-96  to  6-30-97. 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE 

14  July  1997 


4.  TITLE  AND  SUBTITLE  Feasibility  Assessment  and 
Optimization  Study  of  Smart  Actuation  Systems  for  Enhanced 
Aircraft  Maneuver  Performance 


6.  AUTHOR(S} 

K.  Appa,  J.  Ausman,  and  N.S.  Khot 


5.  FUNDING  NUMBERS 


F33615-95-D-3215 

622401 

2401 

55 


WU  00 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AODRESS(ES) 

Northrop  Grumman  Corporation 
One  Hornet  Way 

El  Segundo,  California  90245-2804 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

-  FUGHT  DYNAMICS  DIRECTORATE 
WRIGHT  LABORATORY 
AIR  FORCE  MATERIEL  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433-  7562 
POC:  N.  S.  Khot,  937-255-6992 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  Unlimited 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 

WL-TR-97-3083 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words) 

This  report  presents  an  anlytical  method  developed  to  assess  the  power  requirement 
and  placement  of  the  solid  state  actuators  for  enhanced  maneuver  performance  of 
tactical  fighters.  Optimal  control  methodology  is  employed  to  compute  the  input 
stimuli  such  that  the  total  power  requirement  is  minimum  with  reduced  peak  stress 
levels.  The  study  suggests  that  solid  state  actuators  require  some  form  of 
mechanical  magnification  .of  displacement  for  reduced  control  actuation  power. 


14.  SUBJECT  TERMS 

Smart  Actuation  System,  Feasibility  Study,  Optimal  Control, 
Aircraft  Maneuver 


15.  NUMBER  OF  PAGES 

26 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  I  18.  SECURITY  CLASSIFICATION 
OF  REPORT  I  OF  THIS  PAGE 


Unclassified 


NSN  7540-01-280-5500 


Unclassified 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 


20.  LIMITATION  OF  ABSTRACT 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  Z39-18 


TABLE  OF  CONTENTS 


1  INTRODUCTION  1 

1.1  SUMMARY .  1 

1.2  TECHNICAL  BACKGROUND .  1 

2  ANALYTICAL  MODELING  OF  SOLID  STATE  ACTUATION  SYSTEM  7 

2.1  OVERVIEW .  7 

2.2  STEADY  ROLL  MANUEVER .  7 

2.3  OPTIMAL  CONTROL  DESIGN .  9 

2.4  HAMILTON  EQUATIONS  OF  MOTION .  10 

2.5  SOLID  STATE  ACTUATOR  ELEMENT  FORMULATION .  11 

2.6  SOLID  STATE  ACTUATOR  ENERGY .  12 

3  IMPLEMENTATION  OF  THE  DAMAGE  TOLERANCE  MODULE  13 

3.1  OUTLINE  OF  THE  ALGORITHM .  13 

3.2  USER’S  GUIDE .  16 

3.3  PROGRAMMER’S  GUIDE .  19 

4  DISCUSSION  OF  RESULTS  25 

4.1  TESTS  MODEL .  25 

4.2  TEST  CASE:  STEADY  ROLL  MANEUVER .  26 

5  CONCLUSIONS  AND  RECOMMENDATIONS  31 

REFERENCES  32 


Preceding  Page  Blank 


iii 


PREFACE 


This  report  describes  the  work  performed  by  advanced  structural  methods  department  of 
Northrop  Grumman  Corporation,  Military  Aircraft  Systems  Division,  for  the  Flight 
Dynamics  Directorate,  Wright  Laboratory,  Air  Force  Material  Command,  Wright 
Patterson  Air  Force  Base,  Ohio,  under  contract  F33615-95-D-32 15, Project  No.  0003, 
“Feasibility  Assessment  and  Optimization  Study  of  Smart  Actuation  Systems  for 
Enhanced  Aircraft  Maneuver  Performance.” 

The  final  report  describes  the  theoretical  development  of  the  optimal  control 
methodology  applied  to  determine  solid  state  actuator  placement  and  forces  required  to 
perform  steady  rolling  maneuver  of  flexible  aircraft  at  desired  roll  rate.  This  algorithm  is 
used  employing  ASTROS  generated  data  base.  A  few  test  cases  showing  actuator  power 
requirement  are  demonstrated. 

Dr.  Narendra  S.  Khot,  Structures  Division,  Flight  Dynamics  Directorate,  Wright 
Laboratory,  initiated  the  solid  state  actuator  feasibility  study  program  and  also  served  as 
the  technical  monitor.  Dr.  Kari  Appa  was  the  principal  investigator.  Mr.  John  Ausman 
was  responsible  for  implementing  the  optimal  control  algorithm  in  conjunction  with  the 
ASTROS  module  to  generate  the  required  data. 


1.  INTRODUCTION 


1.1  Summary 

The  performance  characteristics  of  aircraft  largely  depend  on  the  quality  and 
distribution  of  air  flow  on  the  lifting  surfaces.  By  nature,  birds  are  able  to 
configure  their  wings  in  such  a  manner  that  the  air  flow  quality  and  flying 
efficiency  are  at  optimum  conditions.  Duplication  of  such  air  flow 
characteristics  on  manufactured  flying  vehicles  has  been  a  desired  objective  of 
airplane  designers,  beginning  with  the  Wright  brothers.  To  simulate  bird-like 
flying  characteristics,  the  lifting  surfaces  must  be  able  to  deform  smoothly  at 
appropriate  locations.  In  the  mid  1980s,  the  Air  Force  sponsored  a  mission 
adaptive  wing  (MAW)  project  to  study  aerodynamic  and  maneuver 
performance  characteristics  of  tactical  aircraft  (Refs.  1,  2  and  3).  An  F-111 
aircraft  was  selected  and  fitted  with  hydraulic  actuators  to  deform  the  wing. 
This  aircraft  was  test  flown  in  several  mission  performances.  The  test  results 
showed  overwhelming  aerodynamic  performance  benefits  and  agility 
characteristics.  However,  the  actuation  system  was  heavy  and  expensive  to 
operate,  so  practical  implementation  of  this  concept  could  not  be  realized  at 
that  time. 

Recent  analytical  and  wind  tunnel  studies  sponsored  by  ARPA  and  the  Air 
Force  (Ref.  4)  show  how  smoothly  contoured  control  surfaces  promote 
incremental  growth  in  suction  pressure  near  the  leading  edge.  This  has  a 
beneficial  effect  on  control  surface  effectiveness,  leading  to  enhanced  aircraft 
maneuver  performance.  A  few  test  cases  taken  from  this  study  are  presented 
in  the  next  section  to  emphasize  the  need  for  active  control  of  lifting  surface 
camber  to  enhance  aircraft  performance.  This  report,  in  subsequent  sections, 
describes  an  analytical  approach  which  can  be  used  to  determine  optimum 
wing  camber  and  also  to  command  the  solid  state  actuators  to  deform  the 
desired  lifting  surface  to  match  this  optimal  shape. 


1.2  Technical  Background 

In  conventional  aircraft,  the  leading  and  trailing  edge  control  surfaces  are 
used  as  aerodynamic  effectors  to  generate  desired  lift  distributions  on  wings 
and  control  surfaces.  The  leading  edges  are  generally  used  to  minimize  flow 
separation  at  moderately  high  angle  of  attack  flight  maneuver  cases,  while  the 
trailing  edges  are  used  to  obtain  desired  pitching  and  rolling  moments. 

Two  factors  limit  the  effectiveness  of  the  trailing  edge  control  surfaces.  First, 
the  control  effectiveness  decreases  with  increasing  dynamic  pressure  due  to 
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aeroelastic  effects  (i.e.  adverse  twisting  of  the  wing).  The  second  problem  is 
that  there  exists  a  massive  flow  separation  along  the  hinge  line,  resulting  in 
reduced  aerodynamic  loading  on  the  wing  as  well  as  on  the  control  surface. 
Figure  1-la  shows  flow  separation  results  from  a  rapid  change  in  control 
surface  slope  at  the  hinge  line.  This  problem  can  be  alleviated  by  the  use  of 
hingeless  contoured  trailing  edges  as  shown  in  Figure  1-lb  and  Figure  1-lc. 


POINT  OF  FLOW 


POINT  OF  FLOW 
SEPARATION 


(b)  Smart  Contoured  Hinaeiess  Flao 


FLOW  SEPARATION  POIN  T  MOVED 
FURTHER  DOWN  USING  SMART 
ACTUATOR  TECHNOLOGY 


(c)  SMART  CONTOURED  HINGELESS  FLAP  WITH  FLOW 

SEPARATION  SUPRESSION  DEVICES  P96-cm/oi 


Figure  1-1.  Flow  Characteristics  Over  Hinged  and  Hingeless  Contoured 

Control  Surfaces 


A  flow  separation  suppression  device  can  be  used  to  move  the  point  of 
separation  towards  the  trailing  edge.  Use  of  such  devices  will  be  investigated 
at  a  later  time.  To  demonstrate  the  difference  in  pressure  distributions 
between  hinged  flap  and  smoothly  contoured  flaps,  a  few  studies  using  a 
computational  fluid  dynamics  (CFD)  method  were  conducted.  Subsequently, 
wind  tunnel  tests  were  also  conducted  at  NASA  Langley  in  their  Transonic 
Dynamic  Tunnel  (TDT)  under  the  DARPA/WL  "Smart  Materials  and 
Structures  -  Smart  Wing"  contract  (F33615-93-C-3202). 

The  CFD  solutions  are  shown  in  Figure  1-2,  in  which  the  solid  curve 
represents  the  pressure  distribution  computed  for  contoured  trailing  edge 
surface,  while  the  dotted  curve  denotes  the  data  obtained  using  a  hinged  flap. 
The  hinged  flap  data  depict  flow  separation  at  the  trailing  edge,  whereas  the 
contoured  flap  shows  large  suction  pressure  on  the  flap,  as  well  as  on  most  of 
the  upper  surface.  This  type  of  pressure  distribution  has  some  beneficial 
effects  on  aeroelastic  stability.  Since  the  elastic  axis  lies  downstream  of  the 
section  aerodynamic  center,  the  increased  load  near  the  leading  edge  twists 
the  wing  upwards,  resulting  in  an  increased  angle  of  attack  relative  to  that 
observed  in  the  case  of  a  conventional  trailing  edge  flap.  Thus,  the 
effectiveness  of  the  trailing  edge  control  surface  increases  with  increased 
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dynamic  pressure,  and  consequently,  the  roll  reversal  speed  increases.  This  is 
a  significant  contrast  to  the  case  of  traditional  control  surfaces  where  the 
agility  of  the  aircraft  is  reduced  with  increasing  dynamic  pressure. 


Figure  1-2.  Comparison  of  Pressure  Distribution  Between  Hinged  and 
Hingeless  Contoured  Control  Surfaces 


In  the  case  of  wind  tunnel  tests,  two  models  were  constructed.  One  model  had 
conventional  hinged  trailing  edge  flaps  and  ailerons,  while  the  other  (also 
known  as  the  Smart  Wing)  had  deformable  control  surfaces  made  of  shape 
memory  alloys  (SMA).  Typical  wind  tunnel  results  are  presented  in  Figure  1-3 
through  Figure  1-5.  The  trends  are  similar  to  those  observed  in  the  CFD 
solutions.  However,  due  to  lack  of  instrumentation,  the  pressure  loop  for  the 
hinged  control  surface  was  not  observed  in  the  wind  turmel  data. 
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Figure  1-3.  Pressure  Coefficient  (Cp )  Comparison,  at  36%  Span 
Q  =  60  psf,  a  =  8°,  Flap  =  10°,  Aileron  =  0°(  Run  100  vs.  80) 


Figure  1-4.  Pressure  Coefficient  (Cp )  Comparison,  at  50%  Span 
Q  =  60  psf,  a  =  8°,  Flap  =  10°,  Aileron  =  0°(  Run  100  vs.  80) 


‘Smart  Wmg 


smart  Wing 


Conventional  Wing  Conventional  Wingj 


Figure  1-5.  Pressure  Coefficient  (Cp )  Comparison,  at  80%  Span 
Q  =  120  psf,  a  =  8°,  Flap  =  0°,  Aileron  =  0°(  Run  109  vs.  20) 

Aerodynamic  coefficients  (both  lift  and  rolling  moment)  due  to  aileron 
deflection  were  measured  in  the  wind  tunnel  and  were  used  to  compute  the 
roll  reversal  speed.  These  coefficients  are  shown  in  Figure  1-6.  The  roll 
reversal  speed  is  the  velocity  at  which  the  aileron  effectiveness  is  equal  to 
zero.  As  expected,  the  hingeless  control  surface  yields  a  higher  roll  rate  and  a 
higher  control  reversal  speed  compared  to  the  hinged  control  surface. 


0.30 


X=Chord*SQRT(QeCLA/GJ) 

Figure  1-6.  Aileron  effectiveness  Vs.  Nondimensional  Dynamic  Pressure 


Thus,  hingeless  control  surfaces  offer  many  improvements  over  traditional 
control  surfaces.  Chief  among  these  improvements  is  the  improved 
aerodynamic  performance.  Hence,  there  is  a  compelling  reason  to  investigate 
various  avenues,  in  the  light  of  smart  structures  technology,  to  develop 
feasible  mechanisms  to  control  the  camber  of  lifting  surfaces,  as  is  desired  in 
any  combat  mission. 

Today,  with  the  advent  of  new  materials  technology,  it  is  possible  to  design 
smoothly  deforming  lifting  surfaces  using  composite  materials.  Also,  solid 
state  actuators,  which  can  output  large  forces  at  rates  used  in  modern  flight 
control  algorithms,  are  also  being  developed.  These  actuators  can  also  be  built 
in  relatively  small  sizes  and  with  light  weight.  Since  these  actuators  are  small 
and  light,  large  number  of  these  actuators  can  be  used  on  the  lifting  surface  so 
that  any  desired  lifting  surface  deformation  shape  can  be  commanded  for  any 
given  mission  flight  condition. 

To  operate  these  solid  state  actuators  with  minimum  power,  there  is  a  need  to 
determine  appropriate  actuator  locations  and  power  ratings  of  the  individual 
actuators.  This  report  discusses  a  mathematical  approach  based  on  optimal 
control  theory.  Aircraft  performance  goals,  such  as  pitch,  roll  and  yaw  rates, 
are  used  as  the  target  quantities  (constraints),  while  actuator  power  rating  is 
taken  as  the  objective  function  of  the  design  problem.  Detailed  discussions  of 
the  synergetic  design  methodology,  including  the  balancing  of  the  aircraft,  are 
presented  in  References  5  through  8.  A  brief  summary  related  to  the  solid 
state  actuator  is  presented  next. 
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2.  ANALYTICAL  MODELING  OF  SOLID 
STATE  ACTUATION  SYSTEM 


2.1  Overview 

The  main  objective  of  this  study  is  to  develop  an  analytical  model  which 
helps  to  select  the  actuator  locations  and  power  required  for  each  actuator,  so 
that  desired  aircraft  angular  rates  (pitch,  roll  and  yaw)  can  be  achieved 
without  weight  penalty.  A  building  block  approach  is  used  to  develop  the 
algorithm  so  that  more  complex  maneuver  performance  requirements  and 
structural  design  requirements  can  be  added  later. 


2.2  Steady  Roll  Maneuver 

To  develop  the  mathematical  basis  of  the  solid  state  actuation  system,  a  steady 
roll  maneuver  case  is  considered  first.  The  equilibrium  equations  for  this 
problem  can  be  written  as: 

K-r  +  Q-T’^-A-a  +  F-u  =  0  (2-1) 


where: 

Qf  =  Angle  of  attack  at  aerodynamic  panels 

A  =  Aerodynamic  Influence  Coefficients  (AIC)  matrix  with 

respect  to  aerodynamic  panels 

F  =  Nodal  forces  generated  from  the  actuator  elements 

K  =  Structural  stiffness  matrix  in  structural  degrees  of  freedom 
(DOF) 

r  =  Displacement  vector 

T  =  Transformation  matrix  from  structural  DOF  to  aerodynamic 
DOF 

u  =  Vector  of  actuator  stimuli  (input) 

The  displacement  vector,  r,  can  be  expressed  as  a  linear  combination  of  rigid 
body  modes  and  vibration  modes.  Thus,  we  have: 

r  =  [i//,i/rj|jj^|  =  [Vr]77  (2-2) 
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where: 


\j/j.  ~  Vector  of  rigid  body  modes  in  rolling  motion  ((})) 

=  Matrix  of  antisymmetric  vibration  modes 

ri  =  Vector  of  generalized  coordinates,  {rir/Tie} 

(j)  =  Roll  angle  (radians) 

The  subscript  'r'  denotes  rigid  body  modes,  while  'e'  denotes  the  elastic 
vibration  modes. 


The  angle  of  attack  at  the  control  points  of  the  aerodynamic  panels  is  given 
by: 


a  = 

(2-3) 

where: 

V 

=  Free  stream  velocity 

T 

=  Transformation  matrix  relating  the 
structural  points  to  aerodynamic  points 

deformation  from 

P 

=  Roll  rate,  ^ 

u 

=  Chordwise  component  of  V 

The  velocity  component  due  to  rj^  is  omitted  because  it  is  small  in 
comparison  to  the  rolling  velocity,  p. 


Equation  (2-1)  can  be  rewritten  in  terms  of  generalized  coordinates,  T): 


kJ  ^  ^ 


(2-4) 


in  which  the  generalized  stiffness  is  given  by: 


K  = 


0  0 
0  L. 


(2-5) 


This  equation  can  be  solved  for  p  and  1)^: 


=  — [K-F^Aj  3m  = 


(2-6) 


8 


Then  the  roll  rate  is  given  by: 


p  =  Bu  (2-7) 

and  the  generalized  coordinates  of  the  elastic  modes  are: 

Be=Cu  (2-8) 

For  a  given  distribution  of  actuator  stimuli  (u),  one  can  calculate  the  roll  rate 
p  from  equation  2-7.  However,  for  a  specified  roll  rate,  (preqrd)  input 
stimuli,  u,  cannot  be  determined  from  equation  2-7.  This  problem  will  be 
solved  using  optimal  control  theory,  discussed  next. 

2.3  Optimal  Control  Design 

This  problem  can  be  solved  by  the  method  of  optimal  control  theory.  Hence, 
the  Hamiltonian  function  can  be  stated  as: 

h  =  ^  s^Qe  +  UiO  +  (Ap  +  B^)  (2-9) 

where: 

A  =  -1  from  equation  (2-7) 

£  =  Vector  of  target  error  function  which  must  go  to  zero  to 

satisfy  all  performance  requirements 

Q  =  Weighting  matrix 

\  -  Vector  of  Lagrangian  coefficients 

In  the  present  problem,  the  target  vector  (error)  is  given  by: 

a.  Constraint  on  roll  rate: 

8(l)  =  (l-p/p,)  (2-10) 

where  p^  is  the  required  roll  rate  in  the  velocity  range  up  to  . 

b.  Constraint  on  strains  in  the  actuator  materials: 

e(l+j)  =  (l-qu/£3)  (2-11) 
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where: 


=  Allowable  strain  in  the  piezo  electric  material 

Cj  =  Strain  per  unit  of  electrical  input  (volt)  for  actuator 
number  j 

The  main  objective  of  this  constraint  is  to  limit  the  power  input  so  that  the 
actuator  is  not  overly  strained  beyond  elastic  limit.  Likewise,  stress  or  strain 
constraints  on  primary  structural  members  may  also  be  imposed. 

The  second  term,  U,  in  equation  (2-9),  is  the  objective  function  representing 
the  control  power.  It  may  also  include  structural  weight  and  aerodynamic 
figures  of  merit,  which  may  be  considered  in  the  subsequent  studies. 

In  the  following  discussion,  the  state  variable  X  denotes  the  roll  rate,  p, 
and  ^  represents  the  control  variable,  u  (or  the  input  stimuli). 

The  control  power  is  the  sum  of  work  done  by  each  actuator: 

U(^)  =  tw,=^fR§  (2-12) 

I  2 

where: 

^  =  vector  of  control  variables  (u) 

R  =  weighting  matrix 


2.4  Hamiltonian  Equations  of  Motion 


Differentiating  the  Hamiltonian  function,  H  (Equation  2-9),  with  respect  to 

X(=p),  ^(=u),  and  X,  and  using  the  principle  of  optimal  control  theory,  we 
obtain  the  following  two-point  boundary  value  problem: 


(2-13) 


together  with  the  incremental  control  (design)  input: 

^  =  (2-14a) 

in  which 

91  =  [R,+(<2,.,c,/eJ^]  is  a  diagonal  matrix  for  i=l,  n  actuators  (2-14b) 
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The  elements  of  the  Hamiltonian  matrix,  H,  in  equation  (2-13)  are  given  by: 


hyy  —  A 

(2-15) 

(2-16) 

hi  =-QKPt*Pt) 

(2-17) 

^22  ~  1 

(2-18) 

Fy= 

for  i=l,  n  actuators 

(2-19) 

II 

1 

to 

(2-20) 

Since  we  are  dealing  with  a  steady  rolling  maneuver,  the  left  hand  side  of 
equation  (2-13)  can  be  set  to  zero,  and  solved  for  p  and  X.  Thus,  we  obtain: 


(2-21) 


The  actuator  input  for  required  roll  rate  (p^)  is  given  by: 


(2-14a) 


2.5  Solid  State  Actuator  Element  Formulation 

The  actuators  can  be  represented  as  rod  elements,  consisting  of  usual  element 
degrees  of  freedom  and  additional  degrees  of  freedom  to  represent  input 
stimuli  either  in  the  form  of  an  electric  field  or  a  magnetic  field. 

Consider  the  6x1  strain  vector  given  by: 

£  =  8e+£em  (2-22) 


where: 


c  =  Strain  due  to  elastic  deformation 

p  =  Strain  due  to  electric  field 

^em 

can  be  defined  as: 

eeni  =  D-u  (2-23) 
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where: 


u  =  3x1  vector  of  electric  field,  e  (volt /meter) 

D  =  6x3  displacement  matrix  (meter/volt) 


Then,  the  strain  energy  is  given  by: 


U  =  ^\\\e^EedV 


(2-24) 


From  the  variational  principal,  the  element  matrices  can  be  written  as: 


Krr  Ku 

K  K 

^ur  ^uu 


_  \ actuator] 
M j  “  [  sensor  J 


(2-25) 


The  first  set  of  equations  represents  the  function  of  the  actuator  in  which  a 
force  vector  is  generated  due  to  external  stimuli,  u.  The  second  set  of 
equations  denote  the  function  of  a  sensor,  wherein  the  elements  put  out 
electrical  signals  due  to  the  deformation  of  the  structure.  The  stiffness  matrix 
K„  is  computed  in  ASTROS,  while  the  load  matrix  is  computed  in  the  new 
module  and  assembled  into  the  actuator  force  matrix,  F,  as  denoted  in 
equation  (2-1). 


2.6  Solid  State  Actuator  Energy 

Consider  a  piezoceramic  actuator.  The  energy  in  the  element  under 
impressed  electric  field  is  given  by: 


W  =  ^eQe'^e{Vol) 


(2-26) 


where: 


=  Permittivity  of  the  material  (Coulomb  /Newton*Meter^) 
e  =  Electric  field 
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3.  COMPUTER  PROGRAM 
IMPLEMENTATION  WITH  ASTROS 


3.1  Outline  of  the  Algorithm 

The  current  implementation  of  the  Smart  Actuation  system  is  a  two  step 
process,  illustrated  in  Figure  3-1. 


Figure  3-1.  Overview  of  ASTROS/Smart  Actuation  System 
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3.1.1  Step  One 


In  the  first  step,  a  special  version  of  ASTROS  is  run.  This  version  of  ASTROS 
will  process  the  data  that  is  required  for  Smart  Actuation  program.  Two 
additional  bulk  data  entries,  sacntl  and  saelist,  are  required  to  define  data 
that  is  specific  to  this  process.  These  bulk  data  entries  are  described  in 
subsequent  sections.  The  following  files  are  generated  by  ASTROS: 


Format 

Description 

aicmat.d 

1000(1X,E12.5) 

Matrix  of  aerodynamic  influence 
coefficients. 

dyxe.d 

1000(1X,E12.5) 

d¥e. 

dx 

Transforms  elastic  modes  to 
aerodynamic  panels. 

etype.d 

A8 

Element  types 

genkaa.d 

1000(1X,E12.5) 

Generalized  stiffness  matrix. 

geomsa.d 

2(I8,1X),3(1PE12.5,1X) 

Steady  aerodynamics  geometry 

description: 

Column  Description 

1  External  aerodynamic  box  ID 

2  Internal  aerodynamic  box  ID 

3  X  location  of  box  centroid  in 
basic  coordinates 

4  Y  location  of  box  centroid  in 
basic  coordinates 

5  Z  location  of  box  centroid  in 
basic  coordinates 

gtkg.d 

1000(1X,E12.5) 

ASTROS  [GTKG]  matrix. 

parameters.d 

FORTRAN 

'Parameter' 

Problem-dependent  variable  values 
which  should  be  included  in  Kari 

Appa's  Smart  Actuation  program. 

phi.d 

1000(1X,E12.5) 

Matrix  of  mode  shapes. 

phigt.d 

1000(1X,E12.5) 

Transpose  of  ASTROS  [PHIG]  matrix. 

psir.d 

1000(1X,E12.5) 

Rigid  body  modes  matrix. 

xyza.d 

I8,14(1PE12.5) 

Actuator  element  descriptions: 
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The  second  step  of  the  process  is  to  run  the  Smart  Actuation  program.  The 
date  sets  that  were  generated  in  Step  One  are  used  as  inputs  to  the  Smart 
Actuation  program. 


3.2  User's  Guide 


3.2.1  Overview 

This  section  provides  updates  to  the  ASTROS  User's  Manual.  The  sections 
that  are  affected  are  noted  where  appropriate. 


3.2.2  MAPOL  Engineering  Modules 

This  section  provides  additional  documentation  to  Section  2.4.2. 1  of  the 
ASTROS  User's  Manual. 


MODULE 

TYPE 

DESCRIPTION 

SMARTACT 

ENGINEERING 

Compute  and  assemble  data  required  for 

Smart  Actuation  calculations. 

3.2.3  Bulk  Data  Descriptions 

This  section  provides  additional  documentation  to  Section  4.7  of  the  ASTROS 
User's  Manual. 


16 


.nput  Data  Entry  SACNTL 


Description:  Defines  control  parameters  to  the  Smart  Actuation  module. 


Format  and  Example: 


123456789  10 


ITRMAX  Maximum  number  of  design  iterations  (Integer  >  0) 

RMODES  Set  identification  number  of  a  modelist  bulk  data  entry  that  is  used 

to  request  the  rigid  body  modes  that  are  to  be  used  for  the  Smart 
Actuation  analysis  (Integer) 

EMODES  Set  identification  number  of  a  modelist  bulk  data  entry  that  is  used 

to  request  the  elastic  modes  that  are  to  be  used  for  the  Smart 
Actuation  analysis  (Integer) 

ERRTOL  Iteration  convergence  criteria  (Real,  Default  =  0.001) 
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Input  Data  Entry  SAELIST 


Description:  Defines  the  list  of  actuation  elements  to  the  Smart  Actuation 

module. 

Format  and  Example: 


123456789  10 


SAELIST 

SID 

ETYPE 

EIDI 

EID2 

EID3 

EID4 

EID5 

EID6 

CONT 

CONT 

EID7 

EID8 

-etc- 

SAELIST 

1001 

ROD 

1001 

1002 

1003 

1004 

Alternate  Form: 

123456789  10 


SAELIST 

SID 

ETYPE 

EIDI 

THRU 

EID2 

Field 


Contents 


SID  Set  identification  number.  (Integer  >  0) 

ETYPE  Character  input  identifying  the  element  type  (See  Note  1).  One  of 

the  following: 

BAR 

QDMEMl 
QUAD4 
ROD 

E I D  i  Element  identification  number  (Integer  >  0  or  blank) 

Remarks: 


SHEAR 

TRIAS 

TRMEM 


I  Currently,  only  the  rod  element  is  supported  by  the  Smart  Actuation  algorithm.  All 
other  element  types  are  included  here  for  future  development. 

2.  SID  is  here  for  future  implementations  of  the  Smart  Actuation  system  which 
may  be  more  general. 

3.  If  the  alternate  form  is  used,  eid2  must  be  greater  than  or  equal  to  eidi. 

4.  Nonexistent  elements  may  be  referenced  and  will  result  in  no  error  message. 

5.  Any  number  of  continuations  is  allowed. 
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3.3  Programmer's  Guide 


3.3.1  Overview 

This  section  provides  updates  to  the  ASTROS  Programmer's  Manual.  The 
sections  that  are  affected  are  noted  where  appropriate. 


3.3.2  Engineering  Application  Modules 

This  section  provides  additional  documentation  to  Section  5  of  the  ASTROS 
Programmer's  Manual. 
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Engineering  Application  Module:  SMARTACT 
Entry  Point:  SMARTACT 


Purpose: 


To  assemble  and  calculate  data  required  for  Smart  Actuation  analysis 
MAFOL  Calling  Sequence: 

CALL  SMARTACT  (  BC,  SACNTL,  SAELIST,  [PHI6(BC)],  LAMBDA, 

BGPDT(BC),  BEAMEST,  QDMMIEST,  QUAD4EST, 
RODEST,  SHEAREST,  TRIA3EST,  TRMEMEST, 
[DYX] ,  [GENKAA] ,  [GENMAA] , 
AAICMAT(MINDEX) ] ,  TRIM,  OGPWG,  CASE, 
GEOMSA,  [GTKG],  MATl,  MODELIST) ; 


BC 

SACNTL 

SAELIST 
[PHIG(BC) ] 

LAMBDA 

BGPDT(BC) 

BEAMEST 

QDMMIEST 

QUAD4EST 

RODEST 

SHEAREST 

TRIA3EST 

TRMEMEST 


Boundary  condition  number(Integer,  Input) 

Relation  of  smart  actuation  control  parameters 
(Input) 

Relation  of  smart  actuation  elements  (Input) 

Matrix  of  global  eigenvectors  from  real  eigenanalysis 
(Input) 

Relation  of  real  eigenvalue  analysis  results  (Input) 

Basic  grid  point  definition  table  (Input) 

Relation  summarizing  the  cbar  element  (input) 

Relation  summarizing  the  cqdmemi  element  (input) 

Relation  summarizing  the  cquad4  element  (input) 

Relation  summarizing  the  conrod  and  crod 
element  (input) 

Relation  summarizing  the  cshear  element  (input) 
Relation  summarizing  the  ctria3  element  (input) 
Relation  summarizing  the  ctrmem  element  (input) 
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[DYX] 

[6ENKAA] 

[6ENMAA] 

[AAICMAT(MINDEX) ] 

TRIM 

OGPWG 

CASE 

GEOMSA 

[GTKG] 

MATl 

MODELIST 


Elastic  modes  splined  to  the  aerodynamic  model 
(input) 

Generalized  stiffness  (Input) 

Generalized  mass  (Input) 

Antisymmetric  aerodynamic  influence  coefficients 
(Input) 

Relation  containing  trim  parameters  (Input) 

Relation  containing  data  from  the  grid  point  weight 
generation  computations  (Input) 

Relation  containing  the  case  parameters  for  each 
analysis  within  each  boundary  condition  (Input) 

Relation  containing  data  on  the  geometric  location 
of  the  aerodynamic  degrees  of  freedom  (Input) 

Interpolation  matrix  relating  the  forces  at  the 
aerodynamic  degrees  of  freedom  to  the  forces  at  the 
global  structural  degree  of  freedom  (Input) 

Relation  containing  material  properties  (Input) 

Relation  containing  lists  of  normal  modes 


Application  Calling  Sequence: 

None 

Method: 


First  the  case  entries  associated  with  saero  and  saero2  subcases  for  the 
current  boundary  condition  are  read  into  memory.  The  number  of  normal 
modes  is  obtained  by  opening  the  [phig]  matrix.  The  lists  of  normal 
modes  are  read  from  the  modelist  entries.  The  control  parameters  from 
the  SACNTL  entries  are  then  read.  Then  the  list  of  smart  actuation 
elements  is  read  from  the  saelist  entries.  The  material  properties  are 
read  from  the  mati  entries.  The  Basic  Grid  Point  Definition  Table  is  read 
from  the  bgpdt  relation.  For  each  smart  actuation  element,  the  data  in  the 
associated  *est  relation  is  read.  The  G-Set  modes  matrix  is  read  and 
transposed.  The  dyx  matrix  is  read  and  reduced  to  the  number  of  elastic 


21 


modes.  The  gtkg  interpolation  matrix  is  read.  The  phi  matrix  and  xyza 
matrices  are  assembled  for  the  smart  actuation  elements,  phi  describes,  for 
each  mode,  the  generalized  behavior  for  each  smart  actuation  element  and 
has  the  following  form: 


NX 


^PHI  = 


element  1 


element  2 


element  NDVMAX 


~^DOF  of  particular  element 


The  XYZA  matrix  describes  the  smart  actuation  element  connectivity.  The 
generalized  mass  and  stiffness  matrices  are  read.  The  aic  matrix  is  read. 
The  free  stream  velocity  is  read  from  the  trim  relation.  The  center  of 
gravity  is  read  from  the  ogpwg  relation.  The  aerodynamic  geometry  is  read 
in  from  the  geomsa  relation,  psir,  the  rigid  body  aerodynamic  modes 
matrix,  is  assembled.  The  data  required  for  the  smart  actuation  loop  is  now 
assembled. 


Design  Requirements: 
None 

Error  Conditions: 


None 
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3.3.3  Database  Entity  Descriptions 

This  section  provides  additional  documentation  to  Section  9  of  the  ASTROS 
Programmer's  Manual. 


Entity: 
Entity  Type: 
Description: 

Relation 

Attributes: 


SACNTL 


Relation 

Contains  control  parameters  for  the  Smart  Actuation 
System. 


NAME 

TYPE/KEY 

DESCRIPTION 

SID 

Integer  >  0 

Set  identification  number 

ITRMAX 

Integer  >  0 

Maximum  number  of  design  iterations 

RHODES 

Integer 

Set  identification  number  of  a  MODELIST  bulk 
data  entry  that  is  used  to  request  the  rigid  body 
modes  that  are  to  be  used  for  the  Smart 

Actuation  analysis 

EMODES 

Integer 

Set  identification  number  of  a  MODELIST  bulk 
data  entry  that  is  used  to  request  the  elastic 
modes  that  are  to  be  used  for  the  Smart 

Actuation  analysis 

ERRTOL 

REAL 

Iteration  convergence  criteria 

Created  By:  Module  IFP 
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SAELIST 


Entity: 

Entity  Type:  Relation 

Description:  Contains  the  list  of  smart  actuation  elements. 

Relation 

Attributes: 


NAME 

TYPE/KEY 

DESCRIPTION 

SID 

Integer  >  0 

Set  identification  number 

ETYPE 

Text  (8) 

Element  Type.  One  of  the  following: 

BAR 

QDMEMl 

QUAD4 

ROD 

SHEAR 

TRIA3 

TRMEM 

EID 

Integer  >  0 

Element  Identification  Number 

Created  By:  Module  IFF 
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4.  DISCUSSION  OF  RESULTS 


4,1  Test  Model 

To  verify  the  accuracy  of  the  algorithm  discussed  in  previous  sections,  a 
simple  wing  planform,  shown  in  Figure  4-1,  was  selected.  Figure  4-2  shows 
the  smoothly  deformable  trailing  edge  control  surface  mechanism  used  in 
this  analysis.  A  pair  of  solid  state  actuators  were  used  to  deform  the  trailing 
edge  control  surface.  The  actuators  either  pull  or  push  the  stringers.  Thus, 
small  amounts  of  bending  moment  are  applied  to  the  top  and  bottom  skin 
surfaces  so  that  the  control  surface  can  curl  up  or  down  according  to  the 
direction  of  the  stimuli.  In  this  model,  ten  actuators  were  used. 


Figure  4-1.  A  low  aspect  ratio  wing  showing  structural  elements 
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Figure  4-2.  Trailing  edge  control  surface  showing  solid  state  actuation 

mechanism 


A  number  of  performance  analyses  were  conducted.  The  results  are  presented 
next. 


4.2  Test  Case:  Steady  Roll  Maneuver 


Altitude 

Sea  level 

Mach  Number 

LO 

o 

II 

Roll  rate 

Pj=  3.0  radians/  sec 

Number  of  actuators 

10  along  wmg  trailing  edge 

Objective  function 

Minimum  total  power  required 

Constraints 

Target  value  of  p^,  and  strain  allowable  in 
the  actuators 

The  deformed  shape  of  the  wing  is  shown  in  Figure  4-3. 
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The  solid  state  actuators  shown  in  Figure  4-2  are  able  to  magnify  their 
movement,  also  known  as  the  mechanical  advantage  (MA),  such  as  in  the 
case  of  inchworm  motors.  This  is  necessary  to  accommodate  large  relative 
displacements  between  the  end  points  or  attachment  points  of  the  actuators. 
With  this  mechanism,  relatively  small  forces  will  be  able  to  provide  necessary 
wing  camber  for  desired  maneuver  performance.  Table  4-1  through  Table  4-3 
show  the  effect  of  MA  on  actuator  force  requirements  and  consequently  the 
total  energy. 


Mechanical 

Advantage 

1  .OOE-i-02 

Net  Energy 

1017.  ft. lbs 

Actuator 

Stimuli 

Stress 

Force 

Energy 

Displacement 

Element 

Volt 

PSI 

Lbs. 

Inch. lbs 

Inch 

1 

2.22E+01 

-3.17E  +  04 

-1 .59E-I-03 

4.00E-I-01 

0.6 

2 

-2.22E  +  01 

3.91  E-h04 

1.96E-I-03 

6.10E-H01  j 

0.28 

3 

1.52E  +  02 

-2.84E-I-05 

-1 .42E-I-04 

3.21  E-i-03i 

0.495 

4 

-1 .1  2E  +  02 

2.17E  +  05 

1 .09E-I-04 

1 .88E  +  03 

0.026 

5 

1.10E  +  02 

-2.68E  +  05 

-1 .34E-1-04 

2.87E  +  03 

-0.431 

6 

-1 .04E  +  02 

2.49E-I-05 

1.25E-(-04 

2.47E-(-03 

-1  .494 

7 

44.1  9 

-1  33566.6 

-6678.3 

71  1 

-1.52 

8 

-36.21 

1  10627.8 

5531  .3 

488 

-2.93 

9 

29.65 

-83604.3 

-4180.2 

279.1 

-2.28 

1  0 

-24.56 

72489.1 

3624.5 

209.8 

-3.74 

1 

Table  4-1.  Influence  of  Mechanical  Advantage  On  Actuator  Energy 

Requirement,  MA  =  100 


Mechanical 

Advantage 

5.00E+02 

Net  Energy  = 

309. 00ft. lbs 

Actuator 

Stimuli 

Stress 

Force 

Energy 

Displacement 

Element 

Volt 

PSI 

Lbs. 

Inch. lbs 

Inch 

1 

5.2 

-1  0629.8 

-531  .5! 

4.5 

0.598 

2 

-3.7 

1  5760 

7881 

9.89 

0.28 

3 

31.17 

-1  3241  8.1 

-6620.9! 

698.2 

0.489 

4 

-21.57 

103784.2 

5189.2 

428.9 

1 .941 

5 

22.8 

-158524.4 

-7926.2 

1002.2 

-0.435 

6 

-19.9 

143892.8 

71  94.64I 

825.7 

-1  .498 

7 

9.6 

-90266 

-45^3.3 

324.94 

-1.514 

8 

-6.4 

73313.2 

3665.7 

214.34 

-2.921 

9 

6.7 

-54929.37 

-2746. 5l 

120.5 

-2.262 

1  0 

-4.096 

46892.9 

2344. 64| 

87.79 

-3.73 

Table  4-2.  Influence  of  Mechanical  Advantage  On  Actuator  Energy 

Requirement,  MA  =  500 
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Mechanical 

Advantage 

1.00E-1-03 

Net  Energy  = 

251.0  ft.lbs 

Actuator 

Stimuli 

Stress 

Force 

Energy 

Displacement 

Element 

Volt 

PSI 

Lbs. 

Inch. lbs 

Inch 

1 

3.1 

-8150.1 

-407.5 

2.65 

0.59 

2 

-1.33 

1291 1.8 

645.6 

6.64 

0.27 

3 

16.03 

-1  13551.6 

-5677.6 

513.4 

0.48 

4 

-10.24 

89644 

4482.2 

319.9 

0.01 1 

5 

11.87 

-144634.7 

-7231 .7 

834.3 

-0.44 

6 

-9.423 

130595.9 

6529.8 

680.2 

-1.5 

7 

5.29 

-84770.3 

-4238.5 

286.6 

-1.51 

8 

-2.7 

68567.6 

3428.4 

187.5 

-2.914 

9 

3.845 

-51337.5 

-2566.9 

105.2 

-2.24 

1  0 

-1.538 

43671.5 

2183.6 

76.14 

-3.71 

Table  4-3.  Influence  of  Mechanical  Advantage  On  Actuator  Energy 

Requirement,  MA  =  1000 


Figure  4-7  shows  the  energy  required  to  perform  a  3.0  radian/second  roll 
versus  the  mechanical  advantage  of  an  actuator.  For  MA  less  than  100,  a  large 
amount  of  input  power  is  required  to  stretch  the  actuator  to  match  the 
relative  displacement  between  the  end  points  of  the  actuator.  Otherwise, 
desired  roll  rate  cannot  be  achieved.  For  actuators  with  MA  greater  400,  the 
required  energy  is  constant  at  250  ft.lbs.  This  means  the  actuator  is  able  to 
freely  accommodate  the  expansion  between  the  attachment  points. 


40000 

38000 

36000 

34000 

32000 

30000 

28000 

26000 

24000 

22000 

20000 

18000 

16000 

14000 

12000 

10000 

8000 

6000 

4000 

2000 

0 

-2000 


■  Energy  Vs  Mech  Adv.  | 


■  .1.Q0. 


-ZOEL, 


JJZCL 


-4.0.0 _ 5ZC _ ezQ _ zxLQ _ aoo _ aoz _ loioo 


Mechanical  Advantage. 


Figure  4-7.  Actuator  Energy  Vs.  Mechanical  Advantage. 
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5.  Conclusions  and  Recommendations 


An  analytical  simulation  algorithm  based  on  optimal  control  theory  has  been 
developed  to  compute  solid  state  actuator  power  and  placement  to  achieve 
improved  aircraft  performance.  This  study  suggests  that  the  solid  state 
actuators  must  be  able  to  travel  large  distances  between  contact  points  so  that 
the  required  power  is  minimized.  The  displacements  between  points  can  be 
on  the  order  of  half  to  one  half  inches.  Such  large  displacements  and  large 
forces  can  achieved  by  inchworm  actuators. 

The  present  study  used  a  single  degree  of  rigid  body  roll  motion.  Albeit,  this 
algorithm  can  be  extended  to  five  degrees  of  freedom.  Furthermore,  modal 
degrees  of  freedom  can  be  replaced  by  user-selected  structural  degrees  of 
freedom. 
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